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Abstract 

Linear stability of solitary waves near transcritical bifurcations is analyzed for the general- 
ized nonlinear Schrodinger equations with arbitrary forms of nonlinearity and external potentials 
in arbitrary spatial dimensions. Bifurcation of linear-stability eigenvalues associated with this 
transcritical bifurcation is analytically calculated. Based on this eigenvalue bifurcation, it is 
shown that both solution branches undergo stability switching at the transcritical bifurcation 
point. In addition, the two solution branches have opposite linear stability. These stability 
properties closely resemble those for transcritical bifurcations in finite-dimensional dynamical 
systems. This resemblance for transcritical bifurcations contrasts those for saddle-node and 
pitchfork bifurcations, where stability properties in the generalized nonlinear Schrodinger equa- 
tions differ significantly from those in finite-dimensional dynamical systems. The analytical 
results are also compared with numerical results, and good agreement is obtained. 

1 Introduction 

The generahzed nonlinear Schrodinger equations considered in this article are a large class of 
Schrodinger-type equations which contain arbitrary forms of nonlinearity and external potentials 
in arbitrary spatial dimensions. This class of equations are physically important since they include 
theoretical models for nonlinear light propagation in refractive-index-modulated optical media [H 
[2] and for atomic interactions in Bose-Einstein condensates loaded in magnetic or optical traps 
[3] as special cases. Given their physical importance, these equations have been heavily studied 
in the physical and mathematical communities. These equations admit a special but important 
class of solutions called solitary waves, which are spatially localized and temporally stationary 
structures of the system. These solitary waves exist for continuous ranges of the propagation 
constant. At special values of the propagation constant and under certain conditions, bifurcations 
of solitary waves can occur. Indeed, various solitary wave bifurcations in these equations have been 
reported. Examples include saddle-node bifurcations (also called fold bifurcations) [2l|4l[5l[6l[7l[8l|9], 
pitchfork bifurcations (sometimes called symmetry-breaking bifurcations) [3 [lOl [HI [12l [131 [H] > and 
transcritical bifurcations [15j. These three types of bifurcations have also been classified in |15j . 
where analytical conditions for their occurrence were derived. 

Stability of solitary waves near these bifurcations is an important issue. In finite-dimensional 
dynamical systems, stability of fixed-point branches near these bifurcations is well known jl6j . 
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However, as was explained in jSJ [T7j, those stability results from finite-dimensional dynamical 
systems cannot be extrapolated to the generalized nonlinear Schrodinger equations. Thus this 
stability in the generalized nonlinear Schrodinger equations has to be studied separately. For 
saddle-node bifurcations of solitary waves, this stability question has been analyzed in [HI [9]. It was 
shown that no stability switching takes place at a saddle-node bifurcation, which is in stark contrast 
with finite-dimensional dynamical systems where stability switching generally takes place [16] . For 
pitchfork bifurcations of solitary waves, this stability has been analyzed in |1H \T2\ [T^ 117] . It was 
shown that this stability possesses novel features which have no counterparts in finite-dimensional 
dynamical systems as well. For instance, the base and bifurcated branches of solitary waves (on 
the same side of a pitchfork bifurcation point) can be both stable or both unstable [14^ 117]. which 
contrasts finite-dimensional dynamical systems where the bifurcated fixed-point branches generally 
have the opposite stability of the base branch [T6]. For transcritical bifurcations of solitary waves, 
the stability question is still open. 

In this paper, we analyze the linear stability of solitary waves near transcritical bifurcations 
in the generalized nonlinear Schrodinger equations. Our strategy is to analytically calculate bi- 
furcations of linear-stability eigenvalues from the origin when the transcritical bifurcation occurs. 
Based on this eigenvalue bifurcation and assuming no other instabilities interfere, linear stability of 
solitary waves near the transcritical bifurcation point is then obtained. We show that both solution 
branches undergo stability switching at the transcritical bifurcation point. In addition, the two 
solution branches have opposite linear stability. These stability properties closely resemble those 
for transcritical bifurcations in finite-dimensional dynamical systems. Thus, among the three major 
bifurcations (i.e., saddle- node, pitchfork and transcritical bifurcations), the transcritical bifurcation 
is the only one where stability properties in the generalized nonlinear Schrodinger equations closely 
resemble those in finite-dimensional dynamical systems. In the end, we also present a numerical 
example which confirms our analytical predictions. 

2 Stability results for transcritical bifurcations of solitary waves 

We consider the generalized nonlinear Schrodinger (GNLS) equations with arbitrary forms of non- 
linearity and external potentials in arbitrary spatial dimensions. These equations can be written 
as 



where = d"^ /dx\ + d'^ jdx\ + • • • + is the Laplacian in the A^-dimensional space x = 

(xi, X2, • • • ,xn), and F{-,-) is a general real- valued function which includes nonlinearity as well as 
external potentials. These GNLS equations include the Gross-Pitaevskii equations in Bose-Einstein 
condensates [3j and nonlinear light-transmission equations in linear potentials and nonlinear lattices 
[H [21 [H] as special cases. Notice that these equations are conservative and Hamiltonian. 

For a large class of nonlinearities and potentials, this equation admits stationary solitary waves 



where u(x) is a real and localized function in the square- integr able functional space which satisfies 
the equation 



iUt + V^U + F{\U\'^,:x.)U = 







(2.3) 
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and /u is a real-valued propagation constant. In these solitary waves, /x is a free parameter, and u{x) 
depends continuously on fi. Under certain conditions, these solitary waves undergo bifurcations 
at special values of ii. Three major types of bifurcations, namely, saddle-node, pitchfork and 
transcritical bifurcations, have been classified in |15j . In addition, linear stability of solitary waves 
near saddle- node and pitchfork bifurcations has been determined in [8l |9l [TTl [12l [HI [T7]. In this 
paper, we determine the linear stability of solitary waves near transcritical bifurcations in the GNLS 
equations ()2.ip . 

A transcritical bifurcation in the GNLS equations (|2.ip is where there are two smooth branches 
of solitary waves ti^(x;/i) which exist on both sides of the bifurcation point // = ^q, and these 
two branches cross each other at = //q- A schematic solution-bifurcation diagram of transcritical 
bifurcations is shown in Fig. 1(a). Analytical conditions for transcritical bifurcations in Eq. (12. ip 
were derived in [1^. To present these conditions, we introduce the linearization operator of Eq. 

Li = V2-^ + 5„[F(n2,x)n], (2.4) 
which is a linear Schrodinger operator. We also introduce the standard inner product of functions, 

/oo 
r(x)5(x)dx, 
-oo 

where the superscript '*' represents complex conjugation. In addition, we define the power of a 
solitary wave u(x; fi) as 

/oo 
u^(x; fi) dx, 
-oo 

and denote the power functions of the two solitary-wave branches as 

P±{^l) = (n=^(x;/i),u^(x;^i)). 

If a bifurcation occurs at ^ = ^uq, by denoting the corresponding solitary wave and the Li 
operator as 

^to(x) = n(x; /io), LiQ = Li|^=^q^ u=uo, 

then LiQ should have a discrete zero eigenvalue. This is a necessary condition for all bifurcations. In 
|15j . the following sufficient conditions for transcritical bifurcations were derived. In this derivation, 
it was assumed implicitly that the function F(ti^,x) is infinitely differentiable with respect to u. 

Theorem 1 Assume that zero is a simple discrete eigenvalue of Lio at /i = /io- Denoting the 
real localized eigenfunction of this zero eigenvalue as ^(x), and denoting 

G{u; x) = F(n2; x)n, G2(x) = dlG\u=u,. (2.5) 

then if 
and 

A = (1 - GsL^o^no, , V')' - {G2,i'^){G2{L^^u^f - 2L^^uo,i^) > 0, 
a transcritical bifurcation occurs at /i = /io- 
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Perturbation series for the two solution branches u^(x; /i) near a transcritical bifurcation point 
were also derived in [15] . It was found that 

n=^(x; fi) = uo(x) + (/i - fj.o)uf{-x.) + 0{(^ - ^o)^}, (2.6) 

where 

= L^^UQ + 6f V', (2.7) 

and 



= - - - ■ (2.8) 



(1 - GsL^^V') ± 

From these perturbation series solutions, power functions P±{f^) near the bifurcation point can be 
calculated. In particular, one finds that 



p;(mo) = p:(^o), 

thus power curves P±{lj) of the two solution branches ii^(x; /x) have the same slope at the bifurcation 
point. Because of this, the two power curves -P±(^) are tangentially touched at a transcritical 
bifurcation. This feature of the power-bifurcation diagram is shown schematically in Fig. 1(b). 
Notice that this power-bifurcation diagram of the transcritical bifurcation looks quite different 
from the solution-bifurcation diagram in Fig. 1(a). 




Figure 1: Schematic diagrams of a transcritical bifurcation: (a) solution-bifurcation diagram 
(plotted is the deviation function u(xo;/x) — mo(xo) versus /x at a representative xq position); (b) 
power-bifurcation diagram. The same color represents the same solution branch. 

The goal of this paper is to analytically determine the linear stability of solitary waves near 
a transcritical bifurcation point. To study this linear stability, we perturb the solitary waves by 
normal modes and obtain the following eigenvalue problem (see [2], pl76) 



where 



C 










V 







w 



Lo = - /X + F(u2 



X 



(2.9) 

(2.10) 
(2.11) 



and Li has been defined in Eq. (j2.4p . In the later text, operator C will be called the linear- 
stability operator. If this linear-stability eigenvalue problem admits eigenvalues A whose real parts 
are positive, then the corresponding normal-mode perturbation exponentially grows, hence the 
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solitary wave u(x) is linearly unstable. Otherwise it is linearly stable. Notice that eigenvalues of 
this linear-stability problem always appear in quadruples (A, —A, A*, —A*) when A is complex, or in 
pairs (A, —A) when A is real or purely imaginary. 

Using the operator Lq, the solitary wave equation (I2.3P can be written as 



Lou = 0. (2.12) 
In particular, when we denote Lq at the bifurcation point as 

-^'00 = -^'o|/i=/^Oi u=uoi 

then 

Loouo = 0, (2.13) 

thus zero is a discrete eigenvalue of Lqq . 

The main result of this paper is the following theorem on linear-stability eigenvalues of solitary 
waves near a transcritical bifurcation point. 

Theorem 2 At a transcritical bifurcation point /i = /ig in the GNLS equation ()2.ip . suppose 
zero is a simple discrete eigenvalue of operators Lqq and Lio, and 

(^,LooV)/0, P4(/io)/0, (2.14) 

where ip is the real discrete eigenfunction of the zero eigenvalue in Lio (see Theorem 1), then a 
single pair of non-zero eigenvalues itA in the linear-stability operator C bifurcate out along the real 
or imaginary axis from the origin when fi ^ fiQ. In addition, the bifurcated eigenvalues A^ on the 
two solution branches n^(x; fi) are given asymptotically by 



(A±)^ ^ ^7(^-/io), Ai^/uo, (2.15) 



where the real constant 7 is 



7= . (2.16) 

("0, -^^00 V') 



A direct consequence of Theorem 2 is the following Theorem 3 which summarizes the qualitative 
linear-stability properties of solitary waves near a transcritical bifurcation point. 

Theorem 3 Suppose at a transcritical bifurcation point fi = fj,o, the solitary wave uq{x.) is 
linearly stable; and when ^ moves away from no, no complex eigenvalues bifurcate out from non- 
zero points on the imaginary axis. Then under the same assumptions of Theorem 2, both solution 
branches n^(x;/i) undergo stability switching at the transcritical bifurcation point. In addition, 
the two solution branches have opposite linear stability. 

Based on this theorem, schematic stability diagrams for a transcritical bifurcation are displayed 
in Fig. 2. The stability behavior in Fig. 2(a) (for solution bifurcation) closely resembles that 
for transcritical bifurcations of fixed points in finite-dimensional dynamical systems [16]. But the 
power-bifurcation diagram (with stability information) in Fig. 2(b) has no counterpart in finite- 
dimensional dynamical systems. 

Note that for positive solitary waves in the GNLS equations (j2.ip . linear-stability eigenvalues 
are all real or purely imaginary (see Theorem 5.2, pl76). In addition, zero is always a simple 
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Figure 2: Schematic diagrams of a transcritical bifurcation with stabihty information indicated: 
(a) solution-stability diagram; (b) power-stability diagram. The same color represents the same 
solution branch, and solid (dashed) lines are stable (unstable). 

eigenvalue of Lqq [20J. In this case, if zero is also a simple eigenvalue of Liq and the solitary wave 
tio(x) at the bifurcation point is linearly stable, then under the generic conditions (j2.14p . Theorem 
3 applies, thus both solution branches m^(x; fi) undergo stability switching at a transcritical bifur- 
cation point, and the two solution branches have opposite linear stability. Such an example will be 
given in Sec. HI 



3 Proofs of stability results 







0. 



(3.1) 



Proof of Theorem 2 First we see from Eqs. ()2.10p and ()2.12p that zero is a discrete eigenvalue 
of C for all /i values. At the bifurcation point /i = /io, we further have LiQip = 0, thus 



Following the same analysis as in [T7| , we can readily show that the algebraic multiplicity of the zero 
eigenvalue in C is four at = /^o and drops to two when < \^ — <C 1, thus a pair of eigenvalues 
bifurcate out from the origin when ^ ^ fiQ. This pair of eigenvalues must bifurcate along the real 
or imaginary axis since eigenvalues of C would appear as quadruples if this bifurcation were not 
along these two axes. Next we calculate this pair of eigenvalues near the bifurcation point /u = /^o 
by perturbation methods. 

The bifurcated eigenmodes on the solution branches n^(x;^) have the following perturbation 
series expansions, 



(x;m) = X](/'"^o)^^'fc( 



A:=0 



u;^(x;^) = X^ifi - /io)^/2 J^(/i - fio)''w^{x), 

k=0 

X^ip) = iA±(/i - ^o)'/' ( 1 + f;(^ - ^o)'A± ] . 



k=l 



We also expand Lq and Li on the solution branches u^(x; fj.) into perturbation series 

Lq = Loo + {fJ'- f^o)Loi + (^ - Mo)^^o2 + • • • ' 
if = Lio + (^ - ^o)-^^n + (^ - Mo)^^f2 + • • • • 



(3.2) 
(3.3) 
(3.4) 



(3.5) 
(3.6) 
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Substituting the above perturbation expansions into the hnear-stabihty eigenvalue problem ()2.9p 
and at various orders of /i — //Q; we get a sequence of linear equations for {vk,Wk)- 

Lio4 = 0> (3-7) 
Lqqw^ = vq, (3.8) 

LiQvf = (Aq ) wo - L^^o, (3.9) 

and so on. 

First we consider the Vq equation (j3.7p . In view of the assumption in Theorem 2, the only 
solution to this equation (after eigenfunction scaling) is 

= i,. (3.10) 

For the inhomogeneous equation (j3.8p . it admits a single homogeneous solution uq due to 
(j2.13p and the assumption in Theorem 2. Since Lqo is self-adjoint and {uo,ip) = for transcritical 
bifurcations (see Theorem 1), the Fredholm condition for Eq. (j3.8p is satisfied, thus this equation 
admits a real localized particular solution Lqq tp, and its general solution is 

-LqoV + Co^o, (3.11) 



where is a constant to be determined from the solvability condition of the wf^ equation later. 

For the vf^ equation (|3.9p . it is solvable if and only if its right hand side is orthogonal to the 
homogeneous solution ip. Utilizing the and solutions derived above and recalling {uo,ip) = 0, 
this orthogonality condition yields the formula for the eigenvalue coefficient as 

Due to notations ()2.5p and the definition (j2.4p for Li, it is easy to see that L^-^ in the expansion 
(13:6]) is 

L± = G2uf - 1, (3.13) 
where is given in Eq. (gZD. Inserting this if-^ into (IHT^ . we find that 



(Ao^)^ = i^i"^- (3.12) 



(A±)2 = ± . (3.14) 

Substituting this formula into (j3.4p . we then obtain the asymptotic expression for the eigenvalues 
(A^)^ as (I2.15P in Theorem 2. This completes the proof of Theorem 2. □ 

Proof of Theorem 3 Under conditions of Theorem 3, when /i moves away from /xq, the only 
instability-inducing eigenvalue bifurcation is from the origin. We have shown in Theorem 2 that 
this zero-eigenvalue bifurcation creates a single pair of eigenvalues whose asymptotic expressions 
are given by Eq. (|2.15p . This formula shows that on the same solution branch (i.e., n"'"(x;^) or 
ti~(x;/i)), if the bifurcated eigenvalues are real (unstable) on one side of /i = fiQ, then they are 
purely imaginary (stable) on the other side of /u = ^o- Thus stability switching occurs at the 
bifurcation point /U = /Uq for both solution branches. This formula also shows that at the same 
fj, value, if the bifurcated eigenvalues are real on one solution branch, then they would be purely 
imaginary on the other solution branch. Thus the two solution branches always have opposite linear 
stability. This completes the proof of Theorem 3. □ 
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4 A numerical example 



An example of transcritical bifurcations in the GNLS equation (|2.1|) has been reported in ^15j. This 
example is 

iUt + U^cc - V{x)U + \U\^U - 0.2\U\^U + Kc\UfU = 0, (4.1) 
where V{x) is an asymmetric double-well potential 

V(x) = -3.5 sech2(x + 1.5) - 3 sech2(x - 1.5), (4.2) 

and Kc ~ 0.01247946. The potential (j4.2p is displayed in Fig. 3(a), and the power diagram of 
this bifurcation is shown in Fig. 3(b). We see that two smooth power branches, namely the upper 
C1-C2 branch and the lower di-d2 branch, tangentially connect at the bifurcation point {no,Po) ^ 
(3.278, 14.36). This tangential connection agrees with the analytical result on the power diagram 
(see Fig. 1(b)). Profiles of solitary waves at the marked ci, C2, di, ^2 locations on this power diagram 
are displayed in Fig. 3(c-d), and their linear-stability spectra are shown in Fig. 3(ci, C2, di, ^2) 
respectively. These spectra indicate that the solitary waves at ci and d2 of the power diagram are 
linearly stable, whereas the other two solitary waves at C2 and di of the power diagram are linearly 
unstable. Thus both the upper C1-C2 branch and the lower di-d2 branch switch instability at the 
bifurcation point, and the C1-C2 and di-d2 branches have opposite linear stability. These numerical 
results confirm the analytical results in Theorem 3 (see also Fig. 2(b)). 
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Figure 3: Stability switching at a transcritical bifurcation in Example (j4.ip . (a) The asymmetric 
double- well potential (|4.2I) : (b) the power diagram; (c) profiles of solitary waves at locations c\ 
(solid) and C2 (dashed) of the upper power curve in (b); (d) profiles of solitary waves at locations d\ 
(dashed) and d2 (solid) of the lower power curve in (b). (ci, C2, di, ^2) stability spectra for solitary 
waves at locations of the same letters on the power diagram (b) . 
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5 Summary and discussion 



In summary, linear stability of solitary waves near transcritical bifurcations was analyzed for the 
GNLS equations (j2.ip with arbitrary forms of nonlinearity and external potentials in arbitrary 
spatial dimensions. It was shown that both solution branches undergo stability switching at the 
transcritical bifurcation point. In addition, the two solution branches have opposite stability. 
Analytical formulae for the unstable eigenvalues were also derived. These analytical stability results 
were compared with a numerical example and good agreement was obtained. 

The above stability properties closely resemble those in finite-dimensional dynamical systems, 
where it is well known that the stability of fixed-point branches near a transcritical bifurcation point 
exhibits the same behaviors as above [T6|. However, this happy resemblance, which we proved in this 
paper, should not be taken for granted. Indeed, it has been shown that for saddle-node and pitchfork 
bifurcations, stability properties in the GNLS equations differ significantly from those in finite- 
dimensional dynamical systems [8l[9l[I7]. For instance, at a saddle-node bifurcation point, there 
is no stability switching in the GNLS equations ()2.ip [8l[9], but any dynamical-system textbook 
would say that such stability switching takes place [16]. Thus it may be more appropriate to view 
this similar stability on transcritical bifurcations in the GNLS equations and finite-dimensional 
dynamical systems as a happy surprise rather than a trivial expectation. 

Another approach to qualitatively study the linear stability of nonlinear waves in Hamiltonian 
systems is the Hamiltonian-Krein index theory [21 ^12211231124^ 125]. In this approach, the number of 
unstable eigenvalues in the linear-stability operator C is related to the number of positive eigenvalues 
in operators Lq and Li under appropriate conditions. Near a transcritical bifurcation point ^ = ^q, 
we can show that the zero eigenvalue in Lio bifurcates out as 



where A.±[^) is the eigenvalue of Li on the n^(x;/i) solution branch. Using this formula, the 
qualitative stability results in Theorem 3 can be reproduced by the index theory (as was done 
in |17j for pitchfork bifurcations). However, this index-theory approach requires more restrictive 
conditions on the spectra of Lq and Li operators [2ll [24l HZ], and it cannot yield quantitative 
linear-stability eigenvalue formula (|2.15p either. 

It should be recognized that transcritical bifurcations in the GNLS equations (j2.ip occur less 
frequently than saddle- node or pitchfork bifurcations. Indeed, in the example (j4.ip . if the seventh- 
power coefficient k is not equal to that special value Kc, this transcritical bifurcation would either 
turn into a pair of saddle-node bifurcations or disappear, depending on whether k is less than Kc 
or greater than Kc- Due to this less frequent occurrence, one might wonder how useful the stability 
results in this paper are. To address this concern, it is helpful to view a transcritical bifurcation 
as the limit when two saddle-node bifurcations coalesce with each other (such as when k — )• in 
the example (14. ip ). In this connection, the stability results obtained in this paper for transcritical 
bifurcations can also be used to help determine stability properties of nearby saddle-node solution 
branches. Thus the stability results in this article can be useful beyond transcritical bifurcations. 




9 



Acknowledgment 

This work is supported in part by the Air Force Office of Scientific Research (Grant USAF 9550- 
12-1-0244) and the National Science Foundation (Grant DMS-0908167). 



References 

[I] Y. S. Kivshar and G. P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals (Academic 
Press, San Diego, 2003). 

[2] J. Yang, Nonlinear Waves in Integrahle and Nonintegrable Systems (SIAM, Philadelphia, 2010). 
[3] L.P. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Oxford University Press, Oxford, 
2003). 

[4] G. Herring, P.G. Kevrekidis, R. Carretero- Gonzalez, B.A. Malomed, D.J. Frantzeskakis, and 

A. R. Bishop, "Trapped bright matter-wave solitons in the presence of localized inhomo- 
geneities", Phys. Lett. A 345, 144-153 (2005). 

[5] T. Kapitula, P. G. Kevrekidis, and Z. Chen, "Three is a crowd: Solitary waves in photorefractive 
media with three potential wells", SIAM J. Appl. Dyn. Syst. 5, 598-633 (2006). 

[6] C. Wang, P.G. Kevrekidis, N. Whitaker, D.J. Frantzeskakis, S. Middelkamp, and P. Schmelcher, 
"Collisionally inhomogcneous Bose-Einstein condensates in double-well potentials" , Physica D 
238, 1362-1371 (2009). 

[7] T. R. Akylas, G. Hwang, and J. Yang, "From nonlocal gap solitary waves to bound states in 
periodic media", Proc. Roy. Soc. A 468, 116-135 (2012). 

[8] J. Yang, "No stability swtching at saddle-node bifurcations of solitary waves in generalized 
nonlinear Schrodinger equations", Phys. Rev. E 85, 037602 (2012). 

[9] J. Yang, "Conditions and stability analysis for saddle-node bifurcations of solitary waves in gen- 
eralized nonlinear Schrodinger equations", chapter in book Spontaneous Symmetry Breaking, 
Self- Trapping, and Josephson Oscillations in Nonlinear Systems, B.A. Malomed, ed. (Springer, 
Berlin, 2012). 

[10] M. Matuszewski, B.A. Malomed, and M. Trippcnbach, "Spontaneous symmetry breaking of 
solitons trapped in a double-channel potential", Phys. Rev. A 75, 063621 (2007). 

[II] E.W. Kirr, P.G. Kevrekidis, E. Shlizerman, and M.I. Weinstein, "Symmetry-breaking bifur- 
cation in nonlinear Schrodinger/Gross-Pitaevskii equations," SIAM J. Math. Anal. 40, 56-604 

(2008). 

[12] A. Sacchetti, "Universal critical power for nonlinear Schrodinger equations with symmetric 

double well potential," Phys. Rev. Lett. 103, 194101 (2009). 
[13] C. Wang, G. Thcocharis, P.G. Kevrekidis, N. Whitaker, K.J.H. Law, D.J. Frantzeskakis, and 

B. A. Malomed, "Two-dimensional paradigm for symmetry breaking: the nonlinear Schrodinger 
equation with a four-weU potential," Phys. Rev. E 80, 046611 (2009). 

[14] E.W. Kirr, P.G. Kevrekidis, and D.E. Pelinovsky, "Symmetry-breaking bifurcation in the 
nonlinear Schrodinger equation with symmetric potentials", Commun. Math. Phys. 308, 795- 
844 (2011). 

[15] J. Yang, "Classification of solitary wave bifurcations in generalized nonlinear Schrodinger equa- 
tions". Stud. Appl. Math. 129, 133-162 (2012). 



10 



J. Guckenheimer and P. Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifurcations 
of Vector Fields (Springer- Verlag, New York 1990). 

J. Yang, "Stability analysis for pitchfork bifurcations of solitary waves in generalized nonlinear 
Schrodinger equations", to appear in Physica D (2012). 

Y. V. Kartashov, B. A. Malomed, and L. Torner, "Solitons in nonlinear lattices". Rev. Mod. 
Phys. 83, 247-306 (2011). 

G. Hwang, T.R. Akylas and J. Yang, "Solitary waves and their linear stability in nonlinear 
lattices", Stud. Appl. Math. 128, 275-299 (2012). 

M. Struwe, Variational Methods: Applications to Nonlinear Partial Differential Equations and 
Hamiltonian Systems (3rd ed.), Springer, Berlin, 2000. [Specifically the paragraph just below 

Theorem B.4 on page 246.] 

T. Kapitula, P. Kevrekidis, and B. Sandstede, "Counting eigenvalues via the Krein signature 

in infinite-dimensional Hamiltonian systems," Physica D 195, 263-282 (2004). 

T. Kapitula, P. Kevrekidis, and B. Sandstede, "Addendum: Counting eigenvalues via the Krein 

signature in infinitedimensional Hamiltonian systems," Physica D 201, 199-201 (2005). 

D.E. Pelinovsky, "Inertia law for spectral stability of solitary waves in coupled nonlinear 

Schrodinger equations," Proc. R. Soc. A 461, 783-812 (2005). 

M. Chugunova and D.E. Pelinovsky, "Count of eigenvalues in the generalized eigenvalue prob- 
lem," J. Math. Phys. 51, 052901 (2010). 

T. Kapitula and K. Promislow, "Stability indices for constrained self-adjoint operators," Proc. 
Amer. Math. Soc. 140, 865-880 (2012). 



11 



